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"Completeness" (i.e. probability conservation) is not usually satisfied in the cumulant 
expansion of the Anderson lattice when a reduced state space is employed for U — > 00. 
To understand this result, the well known "Chain" approximation is first calculated for 
finite U, followed by taking U — -» 00. Completeness is recovered by this procedure, but 
this result hides a serious inconsistency that causes completeness failure in the reduced 
space calculation. Completeness is satisfied and the inconsistency is removed by choosing 
an adequate family of diagrams. The main result of this work is that using a reduced 
space of relevant states is as good as using the whole space. 

1. Introduction 

With the Hubbard operatorsB one can describe complex local states, as well as 
making substantial simplifications in the study of configurations with a very large 
number jjf states by considering only those few states that are relevant to the 
problerrJ3. The cumulant expansion with Hubbard operators has been employed 
to calculate approximate Green's functions (GF), both for the Hubbard modelEltl 
and for the Anderson latticeErH, but probability is not conserved in the Anderson 
lattice when the state with two local electrons at the same site is projected out of 
the four dimensional space of the local states at each site in the infinite repulsion 
case. As this problem does not appear in the Hubbard model when the whole space 
is considered^, it is important to find out whether this is an essential consequence 
of reducing the space. By considering one of the simplest family of diagrams in 
the cumulant expansion, we can explicitly show that the basic problem is caused 
by an inconsistency that is already present in the full space calculation. The main 
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conclusion of the present work is that the results obtained by employing a given 
family of diagrams in the cumulant expansion would not be worse because of the 
space reduction technique employed in conjunction with the Hubbard operators. 
Although the space reduction is not essential in the Anderson lattice, our result 
gives us confidence in the use of that technique when the study of the full space is 
not viable (e.g. with ions like Eu 3+ that has a ground / 7 configuration with 3432 
states) . 

The Anderson lattice is an important tool in the study of strongly correlated 

r stems, and there are several reviews devoted to this and closely related problems 
. The model consist of a lattice with two localized electronic states at each site, 
that are correlated by a strong Coulomb repulsion U , plus a band of uncorrelated 
conduction electrons (c-electrons) that hybridize with the localized electrons (f- 
electrons). 

The cumulant expansion has been employed by several authors to study the 
Ising and the Heisenberg models SHi, while Hubbard extended the method to a 
quantum problem with fermions. He applied the cumulant approach to the model of 
correlated electrons which he had created, and studied with a different technique 0. 
The method he employed consisted in a perturbative expansion around the atomic 
limit, using the hopping as perturbation, and introducing the Hubbard operators 
Xj,ab =\ j, b)(j, a I which transform the local state | a) at site j into the local state 
| b) at the same site. These X-operators make it possible to describe in a simple 
way the projection of a system into the subspace of the more relevant states, but 
they do not satisfy the usual commutation properties of Bose or Fermi operators. 
Nevertheless, it is possible with the use of cumulants, to derive a diagrammatic ex- 
pansion involving unrestricted lattice sums of connected diagrams, which satisfies a 
linked cluster theorem u. Other expansions with cumu lants i have been also applied 
to study fundamental problems in solid state physics IlilOta, and the relation of 
some of those treatments to the one presented here will be briefly commented in 
the final conclusions. The technique employed by HubbardEHi-3 introduces Grassman 
external fields without physical meaning (unlike the magnetic fields of the corre- 
sponding Ising expansions), but this method seems to be the natural extension of 
the usual expansion with Fermi operatoraj. 

In the present work we shall consider a cumulant expansion of the Anderson 
lattice, which is an extensioncrB of the technique employed by Hubbard to study his 
model. The hopping of the conduction electrons as well as the intra-site Coulomb 
repulsion U between the localized (or f-electrons) are included in the unperturbed 
Hamiltonian, and the hybridization is employed as a perturbation. At a given site 
j, the state space of the f-electrons is spanned by four states: the vacuum state 
| j, 0), the two states| j,a) of one f-electron with spin component ah/2 and the 
state | j, d) with two electrons of opposite spin. In the limit of infinite electronic 
repulsion (U — > oo) the state | j, d) is empty, and one can then consider a reduced 
space of states by projecting | j, d) out. In this space, the identity Ij at site j should 
then satisfy the completeness relation: 
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Xjfio + Xj aa + Xj^o — Ij , (1) 

where a = — a, and the three Xj >aa are projectors into the three states | j, 0) and 
| j, a) of the basis. Because of the translational invariance, the statistical averages 
oi Xj aa (denoted by n a ) are independent of j, and from Eq. ([I]) they should satisfy 

no + n a + n„ = 1 , (2) 

a relation called "completeness" in what follows. It has been found that this relation 
is not usually satisfied when | j, d) is projected out of the space of local electrons 
and the n a are calculated with approximate Green's functions (GF)i. A simple ap- 
proximation displaying this behavior is the "Chain Approximation" (CHA) a0@, 
which is the more general cumulant expansion with only second order cumulants. 
This approximation is <&— derivableEj, and has other interesting properties further 
commented in the Appendix. To better understand the problem, we have studied 
the system for finite U and then taking U — > oo, a procedure that will be called 
CHU in the present work. As in the case of the Hubbard model in the full space!, 
the probability is conserved in the Anderson lattice by this procedure. The analysis 
of the method shows that to secure the completeness, the CHU has to employ dif- 
ferent values for the probabilities of the two states | j, a) and | j, a) , which should 
be equal in the paramagnetic state, and this inconsistency is what causes the lack 
of completeness found with the CHA in the reduced space. This is the basic re- 
sult of the present work, that shows that the calculation in the full space satisfies 
completeness through an important inconsistency, that is hidden in the formalism. 

A systematic way of adding a set of diagrams to an arbitrary family, so that 
completeness be satisfied by the n a calculated_ssnth the corresponding GF, has been 
proposed and verified in a number of cases t3. When applied to the CHA, this 
conjecture gives the "Complete Chain Approximation" (CCHA) which satisfies Eq. 

but it is not derivable any more. It seems interesting to repeat the calculation 
for the diagrams of the CCHA employing a finite U and keeping | j, d) in the space 
of local states, and finally taking U — > oo: this procedure shall be called CCHU in 
what follows. One finds that completeness is also satisfied by the CCHU, and that 
the occupation numbers of the two states | j, a) and | j, a) are equal, and coincide 
with the value calculated directly for the CCHA. This last result, as well as the 
details of the calculation, show that the procedure of making simplified models by 
projecting unoccupied states out of the space of local states, can give adequate 
results. 

2. The Anderson Lattice with Finite U 

The Anderson lattice with finite U is described by the Hamiltonian 

" = E ^k,. c iA« + E ^> fib* + T, U fUi'fUi* + 
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£(v(*)exp(ik.Rj) flC ka +H.C) , (3) 

jka 

where C\ a (Cka) is the usual creation (destruction) operator of conduction band 
electrons with wavevector k and spin component ah/ 2, and the n a ( fj a ) are those 
corresponding to the local (f) electrons at site j. To make the connection with the 
Hubbard operators one should substitute 

fja = Xjficr + uXj^ad (4) 

into Eq. (||), where the factor a = ±1 is necessary to obtain the correct phase of 
the states. As the present treatment employs the Grand Canonical Ensemble of 
electrons it is convenient to introduce 

H = H-J^Cl a C± a + ^v a X j , a A , (5) 
I k,<r i a ) 

where Xj, aa is the occupation number operator of state | a > at site j, and v a 
is the number of electrons in that state. The exact and unperturbed averages of 
the operator A are respectively denoted by < A >u and < A >, and it is also 
convenient to introduce 

£ ]-a = E-j.a - \lV a (6) 

and 

because the energies Ej^ a of all ionic states | a > and the energies E-^ of the 
conduction electrons appear only in that form in all the calculations. 

Operators X^ a \, do not satisfy the usual anticommutation (commutation) prop- 
erties of the Fermi (Bose) operators when the two operators are at the same site, 
but in this case it is sufficient to employ the product rules 

Xj,ab Xj. cc [ — 5b, c Xj. a d . (8) 

One has more freedom to define the operation between pairs of operators when they 
are at different sites, and it is convenient to choose properties as close as possible 
to those of the usual Fermi or Bose operators. It is then convenient to say that 
Xj,ab is of the "Fermi type" ("Bose type") when the number of electrons in the two 
states | j, a) and | j, b) differ by an odd (even) number. The best choice is then to 
use anticommutation relations when the two operators at different sites are of the 
"Fermi type" and commutation relations otherwise. 

The Anderson lattice has been studied in the case of U — > oo resorting to a 
cumulant expansion El that employs the last term in H (hybridization) as a pertur- 
bation, and this treatment is an extension of the expansion already employed by 
Hubbard for his model BO. The Hubbard operators Xj ^ were used in this limit to 
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eliminate the doubly occupied states | j, d) from the state space of the f-electrons 
at site j. The projection of H into the reduced space is 

k.a J> 

The method has been generalized to several configurations with a rather arbi- 
trary choice of states 0, and can be used to study the Anderson lattice with finite 
U , described by the Hamiltonian in Eq. (^). The cumulant expansion gives the GF 
of the Matsubara type for imaginary time t 



Xj, a {T) x r . a ,) ) , (10) 

H 



where X, iQ ,(t) = exp (tH) X^ a exp (— rli.) corresponds to the Heisenberg repre- 
sentation, the subindex a = (b,a) represents the transition | a >— >| b > and the 
subindex + indicates that the operators inside the parenthesis are taken in the order 
of increasing t to the left, with a change of sign when the two Fermi- type operators 
have to be exchanged to obtain this ordering. 

The GFs in Eq. ( |To| ) are defined in the interval < r < [3 = 1/T and, because 
of their boundary condition in this variablea, one can associate to them a Fourier 
series with a coefficient for each Matsubara frequency io v = ttv / [3 (where v are the 
even or odd integers, depending on whether the operator X a is of the Bose type 
or of the Fermi type respectively). One can also transform the GF to reciprocal 
space □ and, because of the invariance against time and lattice translations, they 
are proportional to quantities G^J a , (k, lu v ). One can also associate these quantities 
the iio v points of a complex plane in the variable z = Lu + iy and, in the usual way 
one can make the analytic continuation to the upper and lower half-planes of the 
complex frequency z, obtaining a function G aa ,(lt, z) which is minus the Fourier 
transform of the real time GF. By an adequate choice of a and a' it is possible to 
use these functions to calculate the occupation n a of the local state | a) : 

Pba(u) /tM dLO , (11) 

J — oo 

where 

is the Fermi function and 



/ T (z) = {l + cxp(/3z)}- 1 (12) 



iVs k 



p ba {uj) = - lim Im { — } ^ G baiab (k,v + i |e|) } (13) 
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is the spectral density associated to the transition a — (b, a). With the same pb a (oj) 
it is also possible to obtain the occupation rt& of the local state | b) : 

/oo 
p ba (uj) (1 - /tM) duj . (14) 
-oo 

When the state | d > is projected out, the occupation n a can be calculated with 

only Eq. (|ll|) and G 0(7ct0 , while in the whole space one could also use Eq. ( |l4| ) 
— // 

and G adda . These two values of n a are not automatically equal when one uses 
approximate GFs, and this property is the origin of the results obtained in the 
present work. 



2.1. The chain approximation (CHA) 

The family of diagrams which give this approximation is shown in figure [j]. In the 
diagrams of the cumulant expansion @, the circles (vertices) correspond to cumu- 
lants: the filled ones for the f-electrons and the empty ones for the c-electrons. In 
the present expansion the hybridization is the perturbation, and it is represented 
in the diagrams by the lines (edges) joining two vertices. Because of the form of 
the hybridization in Eq. (|^), the edges can only join a c- vertex to an f- vertex. The 
number of edges reaching a vertex gives the order of the cumulant, and only those 
of second order appear in the CHA. All the cumulants of second order coincide 
with the free propagators, and because of Wick's theorem, these are the only non- 
zero cumulants of the c-electrons, while there are cumulants of any even order for 
the f-elpct rons which can be calculated employing a generalized version of Wick's 
theorem WM 



2.1.1. The Anderson lattice for U — > oo in the reduced space 

The f-electron GF has been calculated for the CHA in the reduced space, and 

ff 

using the simpler notation Go<7,<ro(k, z) = Gq^ a Q(k., z) we can write 

GW,o ff (k,z) = - ? Qg- w . (15) 

(z-e 1 (k))(z-e 2 (k)) 

The energies £i(k) and £2(k) are those of the two elementary excitations with wave 
vector k, resulting from the hybridization of a band and a dispersionless band 
of energy £/ = £j a with a reduced hybridization constant y/ D^ a V(k). They are 
given by the two roots of (z — £f){z — ) — |T^(fc)| 2 = 0, where 

D° Qa = (A 00 + X aa ) (16) 

and (A) is the unperturbed average of A. This GF is equal to that of two un- 
corrected bands of energies £/ and gj. hybridized with a parameter \J D% a V(k), 
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showing that the only effect of the correlations in the CHA is to reduce the hy- 
bridization by a factor y/D^. A typical spectral density po a (ui) is shown in figure 
||, and it is easy to show that 

/oc 
dui AkrM = D° 0a , (17) 
-oc 

i.e. no + n<T, which is the total area of the spectral density in the CHA, coincides 
with the unperturbed value of the same quantity. 



2.1.2. The Anderson lattice for U — ► oo in the full space 

For the Anderson lattice with finite U, described by the Hamiltonian in Eq. (^), 
the same family of diagrams shown in figure |l| corresponds to the CHA. One can 
employ the generalized method introduced in ref. □ for a rather arbitrary choice of 
states, and obtain the approximate GFs for finite U. In reciprocal space and for 
the complex frequency z they are 

G °*MM-- (z _ e/)(1 _ _ &) > (18) 

G* dfia {k,z) = — , (19) 

(z — Ef — U)(l — a — o) 

Gto,™%z) = — D ^ ah - , (20) 



a 



Ch d , m {k,z)=— ° d ) ' , (21) 

\z — Sf — U)(l — a — b) 



where 



DUV{k)\ 



2 



(z-£ f ){z-e ka ) 
D° wd \V(k)\ 2 



(22) 



(23) 



(z-e f -U)(z-e k(7 ) 

r% d = (X^ + X dd ) , (24) 

and -Dg CT is given in Eq. (|l6|). However, it must be noted that now we are employing 
< A > to indicate the unperturbed average of A in the full space, i.e. for the 
Hamiltonian of Eq. (||) taken with V(k) = 0. 
The identity in the full space is given by 

A"j,oo + Xj^ aa + Ajctct + Xj^ dd = Ij (25) 
rather than by Eq. (m), and the completeness with finite U then becomes 



no + n a + n w + n d = 1 



(26) 
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Note that completeness is automatically satisfied for the unperturbed case, because 
then no + n a + n-^ + n& = D® a + D^ d = 1, and it is worth analyzing the two limit 
cases U — and U — > oo for the chain approximation. 



The [7 = case Employing 

<< fjai fj a >>= Goafio- + aG^d,Oa + &Gaa,ad + Gad,ad j (27) 

we obtain the exact GF in this limit, and completeness is satisfied. Although 
this result is not surprising, it serves to illustrate the nature of the X operators, 
which contain strong correlations as a result of their own definition. This is clearly 
shown by employing Eqs. jl9|j20| ) to calculate {Xj, Qa Xj^ d ) and {Xj^ d Xj t0a ): these 
averages should be zero because of the rules of product of the X operators (i.e. Eq. 
(||)), but a different value is obtained in the CHA. This result can be interpreted as 
follows: Wick's theorem is satisfied by the Fermion operators fj a and fj a , and all 
the cumulants of order higher than two are zero, so that the CHA gives the exact 
solution for these operators. As Wick's theorem is not satisfied by the X operators, 
all the diagrams with cumulants higher than two are missing from the individual 



GF in Eqs. (18 -|2 1|) , and they are not exact in the CHA. It is by a cancellation 
automatically built in the formalism that Eq. (^7|) gives the exact << / JCT ; f\ a >>, 
although the Gj a , are only approximate. 

The relation Xj^Xj^d = is an example of the "kinematic restriction which 
forbids more than one energy level of a single ion be occupied" E3. The failure 
of that restriction to be satisfied in the average by many approximate GF of the 
Heisenberg and related models is discussed in El A cure of this problem for the 
Heisenberg model with uniaxial anisotropy is proposed in that reference, and we 
shall later return to this property in connection with our problem. 



The U — ► oo case To study this limit one should start with a very large U . In the 
usual case, the band of conduction electrons is close to Ef, and from Eq. ( p3| ) it is 
b —t Q when U — > oo. A fairly symmetrical situation is obtained when the conduction 
band is close to Ef + U, and in this case it follows that a — > when U — » oo. When 
both Ef and ej + U are far from the conduction band and U — > oo, both b — > 
and a — * 0, and the GF coincide with those of the unperturbed problem. In all 
cases the two "cross" GFs, i.e. GWo<r and G Qa ^d 01 Eqs. (|ll| |o|), vanish for any 
value of z when U — » oo, so that the kinematic restriction discussed in reference 123 
is automatically satisfied in this limit of the CHA. 

Let us consider only the case when Ef is close to the conduction electron band. 
The Eq. Jig} ) then becomes 

D[l(z-£ kff ) 
(z-e 1 (k))(z-e 2 (k)) 
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which is just the same Eq. ( |l5| ) obtained in the CHA when the state | j, d) with 
two f electrons is projected out. In the two ways of dealing with | j,d), Eq. ( |l7| ) 
shows that no + n a coincides with the unperturbed value Do a , and the effect of 
the perturbation consist in transferring weights between no and n„. Equation ( pl| ) 
becomes 

DP. 

G^ d (k,z) = - 7 , (29) 

(z - e f - U) 

and the corresponding p&d is a delta function with weight D^ d at ef + U. Applying 
Eq. (0) and Eq. © to pg d we obtain n d = D% d f T (e f + U) and n w = D% d (1 - 
/t(£/ + U)). The relation n-^ + rid = is then valid independently of the value 



of fi, and together with Eq. (17) shows that the completeness in the full space is 
satisfied: n + n a + n-& + n d = D§ CT + Dl d = 1 . Notice that when Ef is close to the 
conduction electron band, Ef + U is very large and nd — D^ d /t(e/ + £/) — » 0, so 
that completeness in the reduced space n + n a + = 1 is also satisfied. Also the 
unperturbed n° — 0, so that from n-& = D% d follows that = n^-, and this value is 
usually different from the n CT calculated by the CHA. This shows that completeness 
is satisfied in the present calculation at the cost of using different values for n a and 
ti^t. As Eqs. (^8|,^9|) are valid for the two values of a, different values of n a are 
obtained when they are calculated with Goa.oa vs. G a d,ad- When the state | j, d) 
is projected out, the Eq. ( |2^ ) cannot be used, and one is then forced to calculate 

with Gow,oWy so that n-& — n a , and completeness can not be satisfied by the 
CHA in the reduced space. The different values of and n a obtained with the 
CHU are shown in figures | and | as a function of T for the following parameters: 
Ej t(T — —0.5, a local hybridization with V = 0.3 and a rectangular conduction band 
centered at = and with a total width 2W = it. The chemical potential in the 
two figures is respectively fi = —0.60 and u = —0.40, and the corresponding dotted 
curves show the crossover, that occurs at low T in the system with V = 0, from the 
maximal occupation n^_~ 0.5 of the local electron to the empty state n° = 0. 

In many worksB0 it is the GF << f jcr ; f] a » rather than the partial GF 
that is calculated, and completeness is satisfied in the full space. As it was shown 
above, this could hide the fact that there is a pair of conflicting values of n a for 
each a, and that different values of n^ and n a are necessary in the paramagnetic 
case to satisfy completeness: the relation Eq.(^ ) would not be satisfied if only one 
of the two possible values were used for both and n a . The correct completeness 
obtained from << fj a ; y. >> masks this basic failure of the CHA, and the 
properties calculated in the CHU would not be more reliable than those obtained 
in the CHA when | j, d) is projected out. 



3. The Complete Chain Approximation (CCHA) 

A conjecture that gives a systematic way of achieving completeness by adding a 
set of diagrams to an arbitrary family was stated in a previous paper @, and it 
was verified in a number of cases that include diagrams with infinite fourth-order 
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cumulants, but a general derivation has not been found. This conjecture gives for 
the CHA the "Complete Chain Approximation" (CCHA) which satisfies Eq. (S), 
and the extra diagrams are shown in figure |^. We shall now apply to this family 
of diagrams the procedure (named CCHU in this case) employed in Section , of 
keeping | j, d) in the space of local states for finite U and then taking U — > oo. As 
completeness in the full space is satisfied by the CHU, one wonders whether adding 
the extra diagrams to this approximation would maintain Eq. (|26|). This equation 
is satisfied in the CHU at the cost of having different occupation numbers for the 
two states | j, a) and | j,a), and one also wonders whether this inconsistency would 
be removed in the CCHU. In this section it is shown that the two questions have an 
affirmative answer, and that the corresponding n a and n-^ coincide with the value 
calculated directly in the CCHA. 

The fact that different average values of the same operator can be obtained from 
different GF thatxorrespond to the same approximation has been well known in 
magnetic models E3, and in that reference they propose that when a given correlation 
function is calculated, the diagrams should be classified according to the number 
of momentum variables summed over in the diagram and be consistently included 
in the calculation according to this classification (see also references l3cj). This 
procedure is closely related to the conjecture stated in but this last gives clear 
rules for the diagrams that should be included in a GF for theAnderson lattice as 
well as for other systems of fermions. The procedure stated in Ej associates a given 
family of diagrams to the calculation of a given average rather than to a GF: when 
their prescription is applied to our problem, we have a GF in the CHA without 
momentum integrations, but it is necessary to integrate once over the wave vectors 
when the number of electrons is calculated, so that further diagrams should be 
included. Our conjecture seems more clear because it applies directly to the GF, 
and has the same effect as the proposal in @. 



3.1. The GF in the CCHU 

The GFs obtained for finite U in the CCHA are rather complicated, and only the 
corrections AG Q . Q '(k, z) that could give a non-zero contribution to the occupation 
numbers when U — > oo are given below. The calculation of these quantities follows 
the same technique employed in the Appendix C of reference @ to derive the function 
S„ (k, uj u ), which corresponds to the present AGocr,o<r(k, w„) when U — ► oo. We 
obtain: 

AGoa.o^k,^) = - |U(fc)| 2 {A KoaM +B K^(uj u ) C a (cj v )+ 

B' K%Jlo„) H Qd {u v ) + CKl a {Lo u )} , (30) 



AG MiM (k,w„) 



= - \V{k)\ 2 {-A K wd (u„) + B Kl d {u v ) C a (uj v )+ 

B' KU"v) HodM - CKU"u)} , (31) 
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and we do not give the expressions for GVd,oo-(fc, z) and Goa-,wd(k, z) because, as in 
the CHA, they do not contribute to n a when U — > oo. The following definitions 
have been used: 

K ab {u v ) = —± , (32) 

iu y + (e a - e b ) 

N ^ {iuj v - ex) {iuj v - e 2 ) 
k 

PH od {u J ,) = J2c<r(^i)Ko d (u; lJ +Lu 1 ) , (34) 

Ul 

and both e\ = £i(k) and e 2 = £2(k) have been given just after Eq. (|l5|). We 
have left out from the constants in Eqs. ( |30| , pT| ) all those terms that would vanish 
when U — > oo (note in particular that (Xdd) = 0): they can then be expressed by 
A = (X Q0 ) (X aa ) [3 h - D° 0a h;B = - (X„) (1 + D° 0a ); B' = D° od = - (X Q0 ) 
and C = — Z?o CT Ji, where 

/ ( = ^C f W<H (with £=1,2) . (35) 

These values of A, B ancLC are equalElto the parameters with the same name 
in Eq. (4.5) of reference while the parameter B' only appears in the present 
treatment. The terms with A, B and C in Eqs. (30]3l|) are of the same type of the 
corresponding ones in reference E3, and the usual analytic continuation from their 
values at all the points ?„ = iuj v to the complex plane q give analytic functions off 
the real axis. 

The sum over u)\ in Eq. (|34|) was obtained by the usual technique, and the 
resulting Ho d (uj v ) has several terms in which fa(£d — i uiu) is present. To obtain an 
analytic function off the real axis by analytic continuation of the values of Ho d (oj v ), 
it is necessary to use the following property of fa(z): 

fa (e d - i u v ) = (1 - exp (Psdjr 1 = -b T {e d ) , (36) 

where br (z) is the Bose equivalent to the fa (z). It then follows that 

fa (ei) ~ fa ie 2 ) 



ex - e 2 

iu u -e d + e f ( fa (ex) + b T (e d ) fa (£2) + b T (e d ) 



£1 - £2 L 1 u v - e d + £1 1 uj u - £ d + £ 2 

7 / \ i U v - Ed + £f , , 

+br (£d) y. — ry. — r , (37) 

[l UOu - £ d + £1) [l UJ V - £ d + £2) 



tBut note that there is a typing mistake in that reference: the expression for B should read 
B = -(D° (1-D°)+ X „). 
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and the analytic continuation of the total GF would also be an analytic function off 
the real axis, making it possible to employ Eqs. ( pd[]l^ ) to calculate the occupation 
numbers. 



3.2. The occupation numbers for the Anderson lattice in the CCHU 

ft is possible to obtain partial results that are fairly independent of the system 
parameters by integrating Eqs. ([ll 14) over li for each value of k, and leaving the 
sum over this variable as a final step. One calculates the occupation number n a (k) 
associated to the wavevector k and the state | a) by substituting the local spectral 



density p ba (u) in Eq. flllj) or Eq. (|14|) by 



p 6a (k,w) 



(38) 



lim /m jGba,ab( k ^ + ikl) 
Defining Z? <r(k) = ng(k) + n CT (k) and /^(k) = n^-(k) + n d (k) we obtain that 

Z? 0ff (k)=D[] CT -|V(k)| 2 A (39) 

and 



so that 



Z) Fd (k)=D0 d + |V(k)| 2 A 



no(k) + n CT (k) + n^(k) + n d (k) = + Dj d = 1 



(40) 



(41) 



The sum of Eq. (41 ) over all the N values of k divided by N shows that completeness 
(i.e. Eq.(|2^)) is satisfied by any set of system parameters in the CCHU, in particular 
by any unperturbed density p°( £ ) °f c-electron states with respect to their energy 
and by any V{k). 
Using Eqs. (p8|,[30|) we find that 



Mk) = D° Qo 



£i/t (ei) - £2./t (e 2 ) 



h (ei) - h (^2) 



E\ — £2 £l — £2 

\V(k)\ 2 {A f T ( £f ) + B g(£ ka ) + Cf3f T (£ f ) (1 - f T (£ f ))} 



(42) 



where 
5( £ kJ = 



Met) 



h (£1) 



h {£2) 



(e/ - £1) (e; - £ 2 ) (ei - e f ) (ei - e 2 ) (e 2 - £/) (e 2 - £1) 
and employing Eq. ( |39| ) one obtains 

n (k) = D 0ct - n CT (k) 
From Eqs. (p^,pl|) follows that rid(k) = 0, and therefore 



(43) 



(44) 



Mk)=D_ d (k)=D° d + |V(k)| 2 A 



(45) 
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From these expressions is not obvious a priori that one would obtain = rig- 
after summing the Eqs. ( [42| , f45| ) over k, and this property was tested by numerical 
calculation with a rectangular p° c {e). The agreement was perfect, indicating that 
one should be able to prove analytically this property, but we have not made any 
attempt in such direction. A different approach is to compare Eqs. ( |39| , [f2| ) with 
e corresponding expressions that were obtained in the calculation of the CCHA 
: it was verified that the expressions for n a (k) and for Doo-(k) are identical and 
therefore no, n a and D{ )(y = tiq + n a are also equal in the two methods. Since there 
is not an independent expression for in the CCHA, the paramagnetic condition 
n— = n a was imposed E3 and completeness then reads as 

D Qa + n a = 1 , (46) 

which was shown @ to be satisfied for several p° c {e) and V(k). Since Do a and n a 
are identical in the two methods, this relation is then satisfied in the CCHU when 
is true in the CCHA. 

Equation (|d|) shows that when U — > oo, completeness (viz. Eq. (p^)) is satisfied 
for any p° c {e) and V(k) in the full space, and employing n,j(k) = (and therefore 
rid = ) there follows that 



D 0a +n w =l . (47) 

Taking the difference of Eq. (|4^) and Eq. (|47j), it follows that = n a in the 
CCHU and that these values coincide with those calculated with the CCHA, at 
least for all the p° c {e) and V(k) for which Eq. (|^) was numerically verified for 
this approximation. This result shows that employing the conjecture proposed in 
reference the results obtained by the CCHA in the reduced space, that had the 
state | j, d) projected out, are compatible with those calculated by the CCHU in the 
full space. The occupation number of the two spin states is unique in the CCHU, 
thus removing the difficulty presented by the CHA. 

The inconsistency that is brought to light by the use of the CHU is the difference 
between the values of and n a obtained in that calculation. The dashed curves 
in figures || and || plot the value of n a calculated with the CHU, which we have 
shown to be equal to the n a obtained in the CHA, and it is set to be equal to 
the in this last approximation because that quantity can not be independently 
obtained in the reduced space. The calculation in the full space makes possible to 
obtain an independent value of in the CHU, which as discussed before is equal 
to the unperturbed n° , that is shown in figures H and | by the dotted curves. The 
different occupation of and n a in the CHU is clearly shown at low T in the two 
figures, and these are the values that have to be employed to satisfy completeness 
(viz. Eq. ( p6| ) in the full space and Eq. (|^) when U — > oo). It is then clear that it 
is not possible to satisfy Eq. (||) in this region when the relation n-^ = n a is forced 
by the CHA. The n a in the CCHU coincides with the n a in the CCHA, and it is 
given by the full line. 
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It is interesting to remark, that the terms with A and C in Eqs. (|30| and [H]) are 
proportional to (Xoo) {X aa ) /T = n° G n® /T. As this contribution is proportional to 
1/T, it could be dominant at very low T when 0.1 < £//T < 10, but it is very small 
outside this interval because either (Xoo) or (X aa ) is then much smaller than T. 
The region inside this interval is the intermediate valence region, and corresponds, 
at low T, to a crossover from the Kondo region to a region of a thermally excited 
local moment. All the results are fairly well behaved when ej is not too small. 
The whole contribution of Eq. (^o| ) is equal to the difference between the dashed 
curve (n c for the CHU) and the full curve (n a for the CCHU) . This contribution is 
proportional to n®, and the two curves coincide for very small T in figure ^| while 
they tend to different values in figure |J because in this case the n° a does not vanish 
for T -► 0. 

The value of Tig- ciS cL function of fi for T = 0.1 is plotted in figure |^ for the same 
approximations used in the previous figures, showing a monotonic increase with (i 
as it could be expected. 

4. Concluding Remarks 

The Anderson lattice is a model that gives a schematic description of many real 
systems with strongly correlated electrons, and it seems important to understand 
better the different methods employed to calculate its properties. In this work we 
consider the cumulant expansion, a perturbative method that has many desirable 
properties, and we employ the Hubbard operators, that make it easy to project the 
space of the system into a subspace of states of interest. The main motivation of 
using the Hubbard operators is that they are very well suited to reduce the space 
of states of the system by eliminating those that are no relevant to the properties 
considered. 

There are several perturbative expansions that use this operators, and here 
we shall briefly compare three of them: mL3o. All these expansions employ the 
Matsubara techniques with imaginary times, and first expand the grand partition 
function or the Green's functions as a series of statistical averages of time ordered 
products of operators that contain successive orders of the interaction Hamiltonian. 
The technique developed by Keiter, Kimball and Grewe and reviewed in reference 
@, expresses these averages in terms of a diagramatic expansion that leads to a 
scheme of the Brillouin-Wigner tyna^rEZI. This method has been already compared 
with our treatment in more detail El , and shall not be further discussed here. 

We shall consider now the method in E_Hl, and show that it employs a type of 
cumulants different from the ones used in our treatment. The initial averages are 
reduced in E3 by a technique already used by GaudinEa to prove the standard Wick 
's theorem, and later applied to averages of spin operators^ and to Hubbard op- 
erators of the Fermi typajO This technique takes advantage of the fact that the 
commutator of two Hubbard operators is a linear combination of the same type of 
operators, making it possible to reduce an average with n Hubbard operators into 
a linear combination of averages with only n — 1 of them. In 113 they proceed with 
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this reduction until only averages of diagonal Hubbard operators (that are there- 
fore of the Bose type) remain, and these averages are then expressed in terms of 
cumulants that can be obtained from an explicitly given generating function. A di- 
agramatic expansion that employs these cumulants at si ngle, sites is then obtained, 
and a renormalization of vertices can be also performedEj'Eil. 

Our expansion, on thejither hand, closely follows the original derivation em- 
ployed for the Ising model3ol2!a, and the initial averages are expressed directly in 
terms of cumulants, employing a general property that was extended to include 
mulants with Hubbard operators of the Fermi type/ cf. Theorem 3.1 in reference I 
A diagramatic expansion different from that in EL3 is then obtained, and the same 
reduction methoduilj is employed to calculate the expressions of the cumulants. 

Although the two methods are related, the one described in El is different from 
ours, and their cumulants correspond to different operators. The cumulants in El 
only contain "diagonal operators" that are of the Bose type, while the Hubbard 
operators that appear in the cumulants of our expansion are of the Fermi type. In 
particular, the cumulants that are of order n in E3 are of order 2n in our treatment. 
As discussed in the Appendix for the CHA, the relations that are related by Englert 
to the classical Ward identities take in our case a rather different physical meaning 
than the one corresponding to the Ising problem, because in that problem the 
momentum operators are linear combination of Hubbard operators of the Bose 
type, while they are of the Fermi type in our case. 

One difficulty in the cumulant expansion we employ, is that the approximate GF 
obtained do not usually satisfy completeness (i.e. conservation of probability). In 
the present work we consider the cumulant expansion of the Anderson lattice, and 
employ the Hubbard operators to project, in the limit of U — ► oo, the unoccupied 
state with two electrons out of the space of local electrons. It has been shown E£l that 
the "chain approximation" (CHA), which corresponds to the most general family 
of diagrams with only second order cumulants, does not satisfy completeness in the 
reduced space. It was useful to analyze the problem by keeping the doubly occupied 
state for finite U within this approximation, and then taking U — ► oo. It was shown 
that completeness is satisfied when this procedure is used, but to obtain this result 
it is necessary to use different values of n a and n^, that are obtained respectively 
employing Eq. ([ll]) and Eq. ( |l4| ) with the two GFs Go ct ,Oct an d G-^d,ad associated 
with the creation or destruction of the same spin a (cf. Eq. (|27])). To exemplify 
this situation employing either of figures || or ||, note that to satisfy completeness 
it is necessary to use the dotted curve for one spin and the dashed curve for the 
opposite spin. When the GF is calculated in the reduced space (without the state 
| d)) there is no G^d,wd, and the paramagnetic condition n a — has to be forced 
on the system. The dashed line has then to be used for the occupation of both spin 
components, and completeness is therefore not satisfied when the dashed curve is 
different from the dotted one. One important consequence of our derivation is that 
although the calculations employing the full space do satisfy completeness, they 
might hide the use of inconsistent values of n a and n^- 
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A conjecture onhow to obtain families of diagrams that satisfy completeness 
has been proposed E3, and it was verified for some families of diagrams that might 
contain any number of fourth order cumulants E3. The "complete chain approxi- 
mation" satisfies completeness and was obtained by adding diagrams to the CHA 
according to the proposed conjecture. It seemed then relevant to analyze the com- 
pleteness problem in the CCHA by considering again the full space for finite U and 
then taking U — > oo. It was found that the n a and obtained when employing 
the two GFs Go^oo- and G-^d^ad satisfy n a = n-^ and are also equal to the n a that 
was calculated with the CCHA, which only employs Goo-.Oo-- We conclude that the 
inconsistency pointed out for the CHA disappears for the CCHA. Moreover, we 
could find the occupation number n a (k) associated to the wavevector k and the 
state | a) for a — 0, a, a, (cf. Section ) because it is possible to calculate those val- 
ues analytically. Employing Go CTi oct we have also found that the n CT (k) and Uo(k), 
calculated in CCHU are identical to those calculated employing the CCHA. 

The interest in the use of the Hubbard operators is that given a model with local 
states, it is a fairly simple procedure to reduce the corresponding space of states 
to a subspace that only contains the few states that are considered more relevant; 
e.g.: to study the Jahn-Teller effect of Go 2 + in MgO, the 120 states of the (3d) 7 
ground configuration of Co 2+ have been projected into the 12 ground states that 
are split by the Coulomb interaction and by the crystal fields BeII. The question 
immediately arises as to what contributions are lost by the use of this method. The 
technique employed in the CCHU, of comparing the results in the reduced space 
with those obtained by calculating the properties in the full space and then making 
the energy of the states that would be removed go to infinity, seems a good test. 
Although our conclusions strictly apply only to the problem considered, the results 
summarized in the previous paragraphs show that one can safely use the reduced 
space when the occupation of the removed states goes to zero. 
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Properties of the Chain Approximation 

Although the main purpose of the present paper is not to discuss the CHA, and 
much less to propose it for the calculation of all the Anderson lattice properties, this 
approximation has some interesting characteristics that seem worth while presenting 
in a rather summary way. Within the context of the cumulant expansion of the 
Anderson lattice model, the CHA was studied by Hewsoncl and its was recently 
analyzed in more detailEO. 

From the point of view of the high density 1/z classification of diagrams 
the CHA is of zeroth order. An interesting feature is that the CHA is a $— derivable 
approximationEl, generated by the same type of skeleton that gives Jj^c molecular 



field approximation in the IsingLJ and in the quantum Heiscnberg modelQfl. Although 
the present treatment of the Anderson lattice closely follows the cumulant expansion 
for the Ising model a, there is an essential difference with the treatment employed 
for the magnetki^ystems, namely that the external fields we employ to generate 
the cumulantstHallil require the introduction of the Grassman fields £, that have 
no physical meaning and should be set to zero at the end of the calculation. The 
operators associated to those fields are not magnetic moments, which do not change 
the number of electrons, but they create or destroy an electron, so that their averages 
have to be zero when all the £ — 0. As a consequence, one does not obtains a self 
consistent relation equivalent to the Wei^molecular field equation, and for the same 



reason many of the relations established0l§E£l for the Ising model change its meaning 
in the present context. Nevertheless, the CHA retains some of the characteristics 
of a mean field approximationEd: the correlation between the electrons, that in the 
limit U = co would forbid the occupancy of a local electron with spin a at any 
given site that already has a local electron with spin —a, is replaced by an average 
reduction of the hybridization constant V 2 by a factor D a = 1 — n^, equal for all 
the sites. 

When the general formalism of the diagramatic expansion is replaced by an ap- 
proximate evaluation, many of the general properties might be lost in the procedure, 
and one could make a list of properties that would be desirable to maintain E3. It 
seems then interesting to analyze the CHA from that perspective, and in particular 
by considering that it is a $— derivable approximation. 

The study of the <&— derivable approximationsilEl in the cumulant expansions 



^en reviewed by WortisQ, and tha corresponding treatment 
3 and for the Anderson latticala have been also recently 



for the Ising model has 
for the Hubbard model 
discussed, so we shall use the concepts and notations presented in those references 
to avoid repetition. The basic idea is to replace the many unrenormalized diagrams 
of the cumulant expansion by fewer skeleton diagrams with renormalized vertices 
M n , and employ a functional $ of the M n that corresponds to the complete family 
of skeletons, to obtain the grand canonical potential O = T InZ, where the Z is 
the grand partition function. In the derivable approximation only a subfamily 
of all the possible skeletons is considered, and three relations that are valid in the 
exact expansion are employed to define completely the approximate methodLl. These 
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three conditions relate the InZ, M n and the "self-fields" S n and correspond 
to the Eqs.(3.1,3.2,3.5) given in reference £3 for the exact expansion. One would 
then expect that as the correlation functions in a <I>— derivable approximation are 
obtained from "the same underlying free energy" , they "will have singular behavior 
at the same points and obey certain self-consistency conditions"!! 

Variational properties for the free energy that are valid for the Ising modell 
and for the Heisenberg mode£3, can be easily extended to the derivable ap- 
proximations of the Anderson lattice by making the appropriate changes to the 
notatiorfij. These properties guarantee the relation (cf. Eq. (A. 28) in and Eq. 
(3.7) in 0) 



= (?W) c ' =M 1 (£;0 (A.l) 



which corresponds to the Eq. (4.17) in i. Although this relation equates the 
magnetization to the renormalized cumulant Mi in the Ising model, it only gives the 
trivial identity = when £ = in the present cumulant expansion. Nevertheless, 
for £ =/= one can employ 



and, following the Appendix A in reference @, obtain the Eq. (|l5|) employed in 
the present work in the CHA for £ = 0, or the Eqs. for the CHU. One 

could also calculate (Y(£) Y(£)\ directly from the diagramatic expansion, and 

the two results should be equal; Englert has shown that in the Ising model this 
equality requires a relation between the vertex and the bond renormalization ( the 
Eq.(5.22) in reference i), which is very similar to the "Ward identities" of quantum 
electrodynamics. These identities have been also discussed by StmchcomboEl for 
the Heisenberg model. It is a simple matter to show that both the CHA and CHU 
satisfy the Ward identity in the form stated by Englert, because the only non zero 
self- fields in the vertex renormalization are the Si, and the only diagrams that 

contribute to the GF (Y(£) Y(£')) are the simple chains in figure ffl, so that 
the renormalized bond is the same as the bare bond. Another way to prove this 
relation is to notice that the GF in the diagramatic calculation!] are equal to those 
obtained from the $— derivable approximation employing the simple skeleton of 
Figure 3a in reference both for the CHA and for the CHU. 

Another interesting aspect of the CHA is the calculation of the IniJ, that would 
make possible to derive all the thermodynamic properties of the system. By the 
very simplicity of the skeleton that gives the CHA and the CHU, the only non zero 
self-fields in this approximation are the S\l3, and these are also zero in the physical 
region when all £ = 0, so that hiZ = InZo, i.e. we obtain the same grand partition 
function Zq of the unperturbed system. All the thermodynamic relations are then 
automatically satisfied, but at the price of not giving any new information. To 
obtain a better In Z one can use the GF in the CHA, and the coupling constant 
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integration method has been employed to this purpose in reference 0. Another 
possible method is to integrate 

where N t is the total number of electrons in the system. After solving a puzzling 
paradox that appears in a parameter region where three solutions are possible for 
a given N t , one can compare the free energies F obtained by these two methods of 
calculation. In the reference mentioned, their values are plotted for the atomic limit 
and compared with the available exact solution, showing that the two methods give 
fairly close values of F that compare reasonably well with the exact curve. 

Another very important property of the approximation is the spectral density 
of the one-electron GF, which is plotted in Figure |2| for the f electrons, in the CHA. 
This approximation gives the gross behaviour of the spectral density, but it misses 
the presence of the Kondo peak, that should appear close to the Fermi surface at 
low T when the system parameters are in the Kondo region. This peak is what gives 
the characteristically large susceptibility and specific heat of the heavy fermions at 
low T, and its absence from the CHA and CHU is to be expected because in this 
approximation there are no spin flips, that are essential for the existence of the 
Kondo effect. 

To improve on the molecular field approximation one should add the ringsBi, 
that are the next order diagrams in the llz expansion. A family with infinite 
fourth order cumulants has been consideredu, and the corresponding spectral den- 
sity shows some differences with that of the CHA, but it does not show a Kondo 
peak either, as can be seen in the figure 5 of reference □. The same problem is 
therefore present in both the CCHA and the CCHU. One should also mention that 
employing a similar expansion, the magnetization of the Hubbard model has been 
calculated by including all the cumulantsEl, but neglecting the spin fluctuations 
in their calculation: this type of technique would not be then able to recover the 
Kondo peak in the Anderson lattice. By a similajLapproximation the exact solution 
of the Falicov Kimball model has been obtained^, and the spin fluctuations of the 
Hubbard model have been approximately included in the same work. 

The absence of the Kondo peak from the approximate GF that have employed is 
still under study, and we have been able to find an approximation closely related to 
the cumulant expansion that shows a structure near the Fermi surface that decreases 
in intensity when T increaseso, and could therefore be identified as the Kondo peak. 

The present summary shows that the CHA and CHU have very special properties 
within the cumulant expansion, and the absence of completeness in the CHA seemed 
particularly intriguing, thus giving the motivation for the present work, that serves 
to clarify this behavior. 

The CCHA is not $— derivable, but is satisfies completeness and carries within 
its structure the elimination of some inconsistencies of the CHA that are only re- 
vealed through the study of the CCHU. It seems clear that many valuable approx- 
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imations would not satisfy all the properties of the exact solution that one would 
like to maintain, but we believe that they are nevertheless worth its study and use, 
provided that their shortcomings are kept in mind. 
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Figure 1: The diagrams for the GF in the chain approximation (CHA). The circles 
(vertices) correspond to cumulants: the filled ones for the f-electrons and the empty 
ones for the c-electrons. The perturbation (hybridization) is represented in the 
diagrams by the lines (edges) joining two vertices, a) The CHA diagrams for the 
f-electrons GF, represented by the filled square to the right, b) Same as a) for the 
c-electrons, represented by an empty square. 



Figure 2: Plot of the spectral density of f-electrons for the CHA as a function of 
u for a rectangular band of total width 2W = it centered at E = 0, with single 
local electron energy Ej a — —0.5 equal to the chemical potential //, and with 
hybridization parameter V — 0.3, all measured in the same energy units. A gap 
appears at /i w Ef = Ej a — fi. 



Figure 3: The occupation number n a is plotted as a function of T for several different 
approximations and for the same parameters used in figure 0, but with chemical 
potential /i = —0.60. The dashed line is n a calculated with the CHU, and coincides 
with the same parameter calculated with the CHA. The dotted line is calculated 
with the CHU, which coincides with the unperturbed value n° a (i.e. for V = 0). The 
full line is n a calculated with the CCHU and coincides with calculated within 
the same approximation and also with the n a obtained with the CCHA. Note that 
n-a and n a are different for the CHU (dotted vs. dashed line). 



Figure 4: Same as in figure || but for /i = —0.40. 



Figure 5: a) The Complete Chain Approximation: empty squares symbolize the 
diagrams of the c-electron GF in the CHA (cf. figure |l]b) . b) The diagrams added 
to the CHA to make it complete. 



Figure 6: The occupation number n a is plotted as a function of fi for T = 0.1 
while all the remaining parameters are the same used in figure ^ All the curves are 
monotonically increasing functions of fi. 
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